Thermoelectric properties of two-dimensional slabs of BagGaigGeao from first 

principles 
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Thermoelectric effects enable the direct conversion between thermal and electrical energy and pro- 
vide an alternative route for power generation and refrigeration. The clathrate BagGaieGego has the 
highest figure of merit ( ZT ~ 1) among other members in the family of type-I inorganic clathrates. 
Enhancement of the thermoelectric properties have been observed in multilayered superlattices, 
quantum wires and in nanostructured materials, either due to the increase in power-factor (S^a) 
or due to the reduction of lattice thermal conductivity (k). Here, we investigate the thermoelec- 
tric properties of two-dimensional slabs with varying thickness of BagGaigGeso using semi-classical 
Boltzmann transport theory with constant scattering approximation. We observe that, there ex- 
ists a delicate balance between the electrical conductivity and the electronic part of the thermal 
conductivity in reduced dimensions and the insights from these results can directly be used to con- 
trol particle size in nanostructuring experiments. The calculated properties are consistent with the 
recent, first measurements on bulk nanostructured samples. 
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I. INTRODUCTION 



Thermoelectric power generators use the tempera- 
ture gradient to drive the charge carriers (electrons and 
holes) from one end of the material to the other, thereby 
creating a potential difference. This provides the means 
of converting the temperature gradient directly into elec- 
tricity without the need for any moving parts. However, 
the inherent low efficiency of the prevailing materials is 
inadequate to compete with the conventional power gen- 
eration and refrigeration. Hence the central issue in ther- 
moelectrics research is to increase the thermoelectric fig- 
ure of merit ZT, a dimensionless quantity defined as ZT 
= {S'^a / k) /T , where S is the Seebeck coefficient (also 
known as thermopower) , a is the carrier conductivity, 
K is the total thermal conductivity (= k'^^ (electron) 4- 
^p'l (phonon)), and T is the absolute temperature. As 
these transport properties depend on interrelated mate- 
rial characteristics, finding thermoelectric materials with 
enhanced ZT values remains a challenging issue. To en- 
able a more widespread usage of thermoelectric technol- 
ogy power generation and heating/cooling applications, 
ZT of at least 2-3 is required. Two different research 
approaches have been proposed for developing the next 
generation of thermoelectric materials: (i) investigating 
new families of advanced bulk thermoelectric materials, 
and (ii) studying low-dimensional or nanostructured ma- 
terial systems. The increase in ZT of nanostructured 
materials could either be from the reduction of k^^ or 
from the enhancement of the power factor, S'^a. Hence, 
a microscopic picture of the transport properties and its 
response to nanostructuring is necessary to fine tune the 
preparation techniques. 

Type-I inorganic clathrates are host-guest structures 
with the guest atoms trapped in the framework of the 
host structure. From a thermoelectric point of view, the 
clathrates are interesting because they are semiconduc- 
tors with adjustable bandgaps. Investigations in the past 




FIG. 1. (Color online) The crystal structure of type-I 
clathrate BagGaieGeso. The large (green) atoms at the cen- 
ter of the cages are the Ba atoms. The framework is made up 
of two kind of cages: a 20-atom smaller pentagonal dodeca- 
hedron and a 24-atom larger tetrakaidecahedron. The brown 
atoms are Ga and the grey atoms are Ge. 



decade have shown that type-I clathrates ^sGaigGeso 
{A = Ba, Sr, Eu) may have the unusual property of 
"phonon glass-electron crystal" proposed by Slack for 
good thermoelectric materialsJ^ Among the known in- 
organic clathrates, BagGaigGcso has the highest figure of 
merit {ZT « 1). Recently, first successful attempts have 
been made to synthesize nanostructures of BagGaigGeso 
using the bottom-up approach.— The nanoparticles con- 
tained mainly thin plates that were indexed according to 
the [100] direction of the clathrate structure type. The 
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measured Seebeck coefficient reaches -145 /iV/K at 648 
K, though the ZT values are quite low (0.02) due to the 
low sample density. Motivated by these measurements, 
in our work, we concentrate on obtaining a microscopic 
picture of the electronic contributions to the thermoelec- 
tric properties for reduced dimensions. In particular, we 
have considered two-dimensional slabs of varying thick- 
ness with the [100] surface termination and calculated the 
thermoelectric properties using the semi-classical Boltz- 
mann transport theory. We observe that, there exists a 
delicate balance between the electrical conductivity and 
the electronic part of the thermal conductivity in reduced 
dimensions and insights from these results can directly 
be used to control particle size in nanostructuring exper- 
iments. 



II. STRUCTURE 

The type-I BagGaigGcao clathratc has a cubic crystal 
structure with the space group Pm3n. The crystal struc- 
ture contains 46 tctrahedrally coordinated host atoms 
(Ga/Gc) which make two small dodecahedron cages and 
six larger tetrakaidecahedral cages (Fig.[T]). The heavy 
guest atoms (Ba) were placed at the center of the cages 
(wyckoff positions 2a and 6d). The Ga atoms were dis- 
tributed in the framework with twelve Ga atoms at 24fc, 
one Ga at 16i and three Ga at 6c sites, resulting in only 
seven short Ga-Ga bonds with a bond length of w 2.48 A. 
Such a distribution of the Ga atoms was done to accom- 
modate the Ga-Ga bond avoidance observed experimen- 
tally. The remainder of the framework positions were 
then filled with Ge atoms. The symmetry setting for the 
calculations was thus lowered to space group 1 to accom- 
modate such a random distribution of Ga/Gc atoms. 



III. CALCULATIONAL DETAILS 

Non-spin-polarizcd DFT total energy and Kohn- 
Sham band-structure calculations were performed ap- 
plying the full-potential local-orbital code (version 
FPLO9.01.35), within the local density approximation 
(LDA) The Perdew and Wang flavor— of the exchange 
correlation potential was chosen for the scalar relativistic 
calculations. The calculations were carefully converged 
with respect to the number of k points. The transport 
properties were calculated using the semiclassical Boltz- 
mann transport thcorjii^— within the constant scattering 
approximation as implemented in BoltzTraP^iS This ap- 
proximation is based on the assumption that the scatter- 
ing time T determining the electrical conductivity does 
not vary strongly with energy on the scale of kT. Addi- 
tionally, no further assumptions are made on the depen- 
dence of T due to strong doping and temperature. This 
method has been successfully applied to many narrow 
band gap materials including clathrates and as well as to 
oxidesi^i^"— One main concern when calculating trans- 



port coefficients, is the underestimation of band gaps 
using the standard DFT functionals. Such an under- 
estimation of band gaps manifests in the reduction of 
thermopower at higher temperatures due to bipolar con- 
duction. To overcome this problem, we concentrate on 
analyzing the trend in the calculated thermopower for 
various concentrations of electron and hole doping rather 
than quantifying them. 



IV. RESULTS 

As mentioned previously, the nanoparticles of 
BagGaigGeao synthesized recently contained mainly 
nano-plates oriented along [100] direction. Therefore, we 
begin by considering a three-layer slab with a [100] sur- 
face termination and calculate the thermoelectric proper- 
ties. Due to the complex nature of the crystal structure 
of the clathrates, there is no obvious cleavage plane to 
choose from when creating the slab geometries. More- 
over, the cages are all interconnected and every choice of 
termination will destruct the cage network. Hence, we 
calculated the thermoelectric properties of many three- 
layer slabs with different [100] surfaces. Comparing the 
values of Seebeck coefficient for the various calculations, 
we have identified the termination which results in in- 
creasing thermopower as a function of temperature. Hav- 
ing identified an advantageous slab geometry, we pro- 
ceeded to increase the thickness of the slabs and evalu- 
ated the thermoelectric properties without changing any 
other structural parameters. Collected in Fig. [2] are the 
Seebeck coefficient, electrical conductivity and the elec- 
tronic part of the thermal conductivity as a function of 
the chemical potential for a three, four and five layer slab 
for various temperatures. The thermoelectric properties 
of the bulk system is also plotted alongside for easy com- 
parison. The bulk system is a Zintl phase and hence a 
semiconductor. This results in purely electron (chem- 
ical potential ji > 0, S is negative) or hole (chemical 
potential /i < 0, S* is positive) conductivity for low car- 
rier concentrations. On the contrary, the slab geometries 
no longer retain the Zintl composition at the surface and 
hence show multiband behavior, witnessed by the change 
in sign of S for a given /x. 

The values of the calculated thermopower are not 
affected by the constant scattering time approximation 
used to calculate the various transport properties, since 
the expression for S is independent of r. This means 
that S is directly dependent on the electronic structure 
of the material. For the three-layer slab, at /i = 1.2 eV, 
the charge carriers are mainly electrons and S is increas- 
ing with temperature and reaches a value of -90 /iV/K 
at 1100 K. Upon increasing the thickness of the slab to 
four layers, the trend in S is similar and a.t jj, = 1.2 eV, 
it reaches -110 /LtV/K at 1100 K. Adding an additional 
layer retains the large values of thermopower. Compar- 
ing this to the bulk, we can conclude that S remains 
unaffected (even slightly enhanced) upon lowering the 
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FIG. 2. (Color online) Temperature dependence of the calculated thermoelectric properties S, a/r, and ft'^'/r as a function of 
the chemical potential for bulk BagGaisGego along with [100] terminated slabs of varying thickness. 



dimensions for BagGaigGcao. On the contrary, the elec- 
trical conductivity of the slabs reduced by a factor of 
two for the slab systems in comparison to the bulk. In- 
terestingly, the electronic part of the thermal conductiv- 
ity k""' is also reduced vifith respect to the bulk for the 
slab geometries. This term goes into the denominator 
of the expression that defines ZT, and a lower k'^^ for 
the slabs means an increase in ZT values with respect 
to the bulk. Assuming that nanostructuring reduces the 
phononic part of the thermal conductivity k^^ by intro- 
ducing grain boundaries, one can obtain improved ther- 
moelectric performance in nanostructured BagGaigGeig 
clathrates. 



V. CONCLUSIONS 

We have calculated the electronic contribution 
to the thermoelectric properties for two-dimensional 



slabs of BagGaigGcao clathrate from semi-classical 
Boltzmann statistics. We discern a dependence of the 
different surface terminations and slab thickness on the 
calculated thermoelectric properties. A delicate balance 
exists between S, a and k^' in reduced dimensions. 
Enhanced values of S is obtained for heavily n-doped 
slabs of BasGaigGeao with respect to the bulk. Both 
CT and k"^' are reduced for the slabs in comparison to 
the bulk, but S'^a/n'^^ of the slabs are comparable to 
the best bulk values. Assuming a reduction of k^'' via 
nanostructuring, one can anticipate improved values of 
ZT for low dimensional BagGaigGcao systems. 
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